Interaction of Josephson and magnetic oscillations in Josephson 
tunnel junctions with a ferromagnetic layer 
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We study the dynamics of Josephson junctions with a thin ferromagnetic layer F [superconductor- 
ferromagnet-insulator-ferromagnet-superconductor (SFIFS) junctions]. In such junctions, the phase 
difference ip of the superconductors and magnetization M in the F layer are two dynamic parameters 
coupled to each other. We derive equations describing the dynamics of these two parameters and 
formulate the conditions of validity. The coupled Josephson plasma waves and oscillations of the 
magnetization M affect the form of the current-voltage {I-V) characteristics in the presence of a 
weak magnetic field (Fiske steps) . We calculate the modified Fiske steps and show that the magnetic 
degree of freedom not only changes the form of the Fiske steps but also the overall view of the I-V 
curve (new peaks related to the magnetic resonance appear). The I-V characteristics are shown 
for different lengths of the junction including those which correspond to the current experimental 
situation. We also calculate the power P absorbed in the system if a microwave radiation with an 
ac in-plane magnetic field is applied (magnetic resonance). The derived formula for the power P 
essentially differs from the one which describes the power absorption in an isolated ferromagnetic 
film. In particular, this formula describes the peaks related to the excitation of standing plasma 
waves as well as the peak associated with the magnetic resonance. 

PACS numbers: 74.20.Rp, 74.50.-l-r, 03.65.Yz 



I. INTRODUCTION 

A great attention in recent years has been paid to the 
study of Josephson junctions (JJ) with a magnetic layer 
(or layers) ^"^. Ahhough the exchange field in the fer- 
romagnetic layer F essentially suppresses the Josephson 
current /j, the interaction of the exchange field and sin- 
glet Cooper pairs results in new, interesting, and nontriv- 
ial effects. For example, the singlet pair wave function 
penetrating from the superconducting leads into the F 
layer due to the proximity effect oscillates in space. In 
case of a uniform F layer, the pair wave function consists 
of two components: one is the singlet component and 
another is the triplet component with zero projection of 
the total spin on the direction of the magnetization vec- 
tor M in the ferromagnet. The condensate wave func- 
tion decays in the ferromagnet on a short distance from 
the superconductor-ferromagnet (SF) interfaces, which 
in the diffusive limit, is of the order £^p = ■y/_D/2i?oxc, 
where D = vpl/S is the diffusion constant and i?exc is 
the exchange energy. Here, vp and I denote the Fermi 
velocity and the electron mean free path, respectively. 
Oscillations of the Cooper pair wave function in space 
lead to a change of sign of the critical Josephson current 
Ijc- This effect was predicted long ago'"''*' but observed 
only recently ^"^''. 

If the magnetization in the F layer is not uniform (for 
example, this occurs in the case of a domain structure 
or multilayered ferromagnet-superconductor (FS) struc- 
tures with noncollinear magnetization directions in the 
F layers), due to the proximity effect a so-called odd- 
frequency triplet component arises"^'"*'^^. In contrast to a 
conventional triplet component that is an odd function 
of momentum and is suppressed by scattering off ordi- 



nary impurities , the odd- frequency triplet component 
is an even function of momentum (in the diffusive case) 
and is not destroyed by scattering off ordinary impurities. 
This component also is not sensitive to the exchange field 
and therefore can penetrate into the ferromagnet over a 
long distance up to ^^v = \J DJ^-kT at temperature T. 
Convincing data in favor of existence of this long-range 
triplet component have been obtained in a number of 
recent experimental works^^"^"^. 

Another interesting effect arises in SFIFS junctions. 
It turns out that at the antiferromagnetic magnetization 
orientation in the F layers, the Josephson critical current 
Ijc is increased^^. Its value may even exceed the critical 
current Ijc in similar JJs without ferromagnetic layers. 
This prediction was also confirmed experimentally^^. 

Alongside with the study of the dc Josephson cur- 
rent in SFS or SF1F2F1S JJs, dynamic properties of 
these junctions and also of tunnel SIFS or SFIFS JJs 
have been investigated both experimentally'^^"^^ and 
theoretically"^""'^^. Here and throughout the paper, S and 
I, respectively, represent a superconducting and insulat- 
ing layers and F1/2 denotes two distinct ferromagnetic 
layers. Interesting dynamic phenomena in JJs with a fer- 
romagnetic layer or a magnetic particle occur when the 
dynamics of the superconducting phase difference (/9(t) 
and the magnetization Af (i) come into play. 

The coupling between these two degrees of freedom 
may be realized in different ways. For example, the 
Josephson current produces a torque acting on magne- 
tization vectors in multilayered SF1F2S junctions. Since 
the Josephson current /j[(/3(t)] is determined by the mu- 
tual orientations of magnetization vectors M1/2, the dy- 
namic behavior of the Josephson current will depend on 
the dynamics of M(tf'^''^'° . Another mechanism of the 



supercurrent action on magnetization was considered by 
Konschelle and Buzdin'^^. They studied dynamics of SFS 
junctions with a non-centrosymmetric ferromagnet. In 
this case, the Josephson current Ij acts directly on the 
magnetization M leading to its precession. In a nonsta- 
tionary case, the interplay between Ij{t) and M{t) leads 
to a complicated behavior of the phase difference (f{t) in 
time. 

In several papers, '^^"'^^ dynamics of SmS 
(superconductor-magnetic impurity-superconductor) 
JJs have been studied, where m stands for a magnetic 
impurity. Interaction between tunneling Cooper pairs 
and the magnetic moment of the impurity not only 
changes the current-phase relation Ij{tp) but also results 
in interesting dynamics of the magnetic moment. 

The most interesting dynamic effects arise in tunnel 
JJs with a ferromagnetic layer (or layers). In this case, 
the interaction between the magnetization in F and the 
Josephson current is realized in the simplest way. As is 
well known, even a weak in-plane magnetic field strongly 
affects the Josephson current /,/(</?). In case of JJs of the 
SIFS or SFIFS type, such a magnetic field is produced 
by the F layer itself. Therefore any perturbations of the 
magnetization vector M change the current Ij{^) and 
in addition the Mcissner currents in the superconducting 
leads change the orientation of the M vector. 

In absence of the F layer, Josephson plasma waves can 
propagate in SIS junctions and their spectrum is^*'^^^: 
cj^ = rtj + k'^Vj, where fij is the Josephson "plasma" 
frequency and vj is the velocity of Swihart waves. On 
the other hand, in the F film, spin waves can be excited 
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the magnetic resonance frequency and Im is a "magnetic" 
length"'^. If i}j < Hm, then these dispersion curves cross 
(usually ^^j -^ Pj = Vj/ft^j), and the interaction between 
magnetization and Josephson currents leads to a coupling 
between Josephson "plasma" and spin waves and to a re- 
pulsion of the corresponding dispersion "terms." The 
coupling between magnetic and superconducting oscilla- 
tions can be observed by studying the I-V characteristics 
of the junction in the presence of a weak external mag- 
netic field. In this case, the so-called Fiske steps arise on 
the I-V curve, but their particular positions and form de- 
pend on parameters characterizing the magnetic system. 
New peaks related to magnetic resonances appear on the 
current-voltage characteristics (CVC). These results have 
been obtained in a short paper by two of us^°. 

In the current paper, we study dynamic phenomena in 
the same systems (SIFS or SFIFS JJs) as in Ref. 30. 
However, we present in more detail the derivation of 
equations describing the dynamics of the coupled mag- 
netic and superconducting systems (see Sec. II). In par- 
ticular, we formulate conditions (frequency range) under 
which these equations are valid. As in Ref. 30, we ana- 
lyze Fiske steps in SFIFS junctions, but the CVC will be 
presented for a wider range of parameters of these junc- 
tions. The CVC will be displayed not only for junctions 
with L = Ij as it was done in Ref. 30, but for junctions 
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FIG. 1. (Color online) Schematic construction of a SFIFS 
junction of the "overlap" geometry. 



longer or shorter than the Josephson length {L < Ij). 
The latter case corresponds to the current experimental 
situation. 

The coupled magneto-plasma modes will also be dis- 
cussed in more detail (see Sec. IV). Finally, in Sec. V, we 
present a formula for the power absorption V in SFIFS 
junctions when a weak ac in-plane magnetic field is ap- 
plied, that is, we study the ferromagnetic resonance in 
the system. This formula drastically differs from the 
known formula for ferromagnetic resonance in an isolated 
F film. In particular, it describes plasma resonances in 
tunnel JJs, which also occur in absence of the F film. The 
frequency dependence of V will be presented for various 
system parameters. In Sec. VI, we discuss the obtained 
results and analyze possibilities to observe the predicted 
effects in experiments. 



II. MODEL AND BASIC EQUATIONS 

We consider a planar SFIFS junction of the "overlap" 
geometry as shown schematically in Fig. 1 (the results ob- 
tained are also applicable to an SIFS junction). Our aim 
is to generalize the equation for the phase difference tp 
between the superconducting layers describing the static 
and dynamic properties of an SIS J J to the case of SFIFS 
JJs. 

This equation reads^^"^^'**'*'^^ 
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fjVlp + sm{p)=rj, (1) 



where flj = (2ejc/C'nfi)^^^ is the Josephson "plasma" 
frequency, 7^ — {R\jC\j)^^, C\j = e/Aird, and R\j are 
the capacitance and resistance of the junction per unit 
area, respectively, d is the thickness of the insulating 
layer, Zj = Vj/fl'^, vj — c^/d/2eXi, is the plasma wave 
propagation velocity (Swihart waves), Xl is the London 
penetration depth, and V_l represents the tangential or 
in-plane gradient with respect to the interfaces in the x-y 
plane. 



We single out the term on the right-hand side of 
Eq. (1), 7/ = j /jci which describes the normaUzed bias 
current through the junction. Although it may depend 
on y, the normalized current rj will be considered as con- 
stant along the y direction. Strictly speaking, this is only 
true for "overlap" j unctions'^ ^'^^ considered here in which 
the system geometry is arranged in such a way that the 
intersection region of superconducting layers is approx- 
imately one-dimensional. However, the form of Eq. (1) 
is most convenient for analysis of CVC for the system 
under consideration and, moreover, neglecting the y de- 
pendence of normalized current 77 does not change qual- 
itatively the final results. The critical current density jc 
is considered as a known quantity. It was calculated in 
Refs. 24, 46-48. 

The resistance i?n depends on the voltage V across the 
junction. This dependence is especially strong in the case 
of tunnel SIS JJs if the voltage V is close to the energy 
gap A. We assume that the characteristic frequencies 
(Oj and JIm) are smaller than A//i. In addition, we are 
interested in the form of the CVC at voltages V close 
to hQ.j/2e, hflM/^e, where i?n and, therefore, 7^ can 
be regarded as constant. Of course, the overall form of 
the CVC will be modified as a direct consequence of the 
voltage-dependent damping coefficient Jr{V). 

We consider planar JJs of the SFIFS, SFIS, or SFS 
type and assume that the layer separating the two su- 
perconductors is characterized by the magnetic suscep- 
tibility x(a;,fc). In particular, this layer may be a mag- 
netic insulator or metallic ferromagnet. The derivation 
of an equation for the phase difference ip in SFIFS junc- 
tions is quite similar to that in the case of tunnel SIS 
junctions'*^''*^'^^'^^. We assume that there is no mag- 
netic field normal to the interfaces in the superconductors 
or, in other words, no Abrikosov vortices pierce the su- 
perconducting films, and the lateral dimensions L^^y are 
much larger than the thickness dp of the F layers and the 
Josephson penetration depth A^. Since the normal com- 
ponent of the magnetic induction B^ is continuous at the 
superconductor-ferromagnet (SF) interfaces, it also van- 
ishes in the ferromagnetic layers and, hence, according to 
B^ = H:,+ AttM^ one has H,_ = -At:M^ in the F films. 
In order to find the relation between the magnetic field 
H in the superconductor (note that in the S layers H 
coincides with the magnetic induction B) and the phase 
difference (/?, we express the tangential component of the 
current density in the S film jj_ = jx^^x + jy^y using the 
vector potential A± (n^ x dA/dz — 'B±) and the tan- 
gential gradient of the phase in the superconductor Vj^x 



as 
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where jqp{iu) — 47riwcr(aj)A|^/c^ is a damping parameter 
describing effects of quasiparticles on the supercurrcnt 
and ^o = hc/2e > is the magnetic flux quantum. The 
parameter "fqp is very small for not very high frequencies 



because the frequency c/Xl is very large. For example, 
taking Al = 5 • lO'^cm we obtain c/Xl = 0.6 • 10i'^s-\ 
which actually allows us to omit the parameter 7^^. 

Writing Eq. (2) we imply a local relation between the 
tangential current density jj^ and the gauge invariant 
quantity in brackets, which is legitimate in the limit 
kX^ ^ 1, where k is the modulus of the in-plane wave 
vector of perturbations. Subtracting the expressions for 
the current density, Eq. (2), written for the right and 
left superconductors from each other we find the change 
of the tangential current density [j_l] — i±{dF/2) — 
jj_{~dp/2) across the junction 



[J±\ = 
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(3) 
where dp — dp in the case of an SIFS or SFS junction 
and dp = 2dp in the case of an SFIFS junction. The 
parameter dp is the thickness of the F film, which is as- 
sumed to be smaller than the London penetration length 
Al, and for any quantity Q, we denote the difference 
^\s(R) ~ ^\s(L') ^y ['31' where S{R) and S{L) are the 
right and left superconductors, respectively. 

The assumption dp <^ Xl allows one to neglect the 
change of A_l along the z direction caused by Meissner 
currents in the F layer and to write the change of the 
vector potential A_l in the form [Aj^] ~ dp{nz x Bj^) 
with Bj^ = 47rMj_ -|-H_l. The field Hj^ is approximately 
the same to the right and to the left from the SF inter- 
faces and does not contribute to the jump of the tangen- 
tial current density [jj^]. The Meissner currents in the F 
layers and, therefore, the variation of Hj^ there are much 
smaller than in the superconductors for the following rea- 
son. The total screening Meissner current /scr in the F 
layer is proportional to Xj^pdpA, where the inverse Lon- 
don penetration depth Xj^p is proportional to the density 
of Cooper pairs, Xj^p ~ nsp, and, thus, is much smaller 
than A£ . The phase difference (p between the two S 
layers has the (gauge-invariant) definition: 



'P=[X] + 
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S{R) 



dz A, 



(4) 



S{L) 



and completely describes the JJ because we choose a 
gauge with A^ ^ and [x] = x{dp/2) - x{~dp/2). 

Equation (3) determines the boundary conditions of 
the London equation in the superconductors. Indeed, 
considering the Maxwell equation at the points z = 
±zsF ~ ±dp/2, 



V x B_L 
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±dp/2 
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±dp/2 



(5) 



where Zsf denotes the coordinate of the right SF inter- 
face, we obtain by successively taking the cross product 



with riz in both sides and subtracting the two equations 
from each other 



dB^ 



dz 



Here, we used the relation 



27r 



dF/2 



^z X [J_l]. 



(6) 
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Finally, with the help of Eq. (10) and taking into account 
that in the S layers B = H, 



9B_L 



dz 
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dz 
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taking into account the symmetry of the SFIFS system. 
Recalling that the magnetic field component B^ normal 
to the interfaces is assumed to be zero in the S layers 
and considering only the z dependence of Bj^, we have 
to solve in the superconductors the equation 



dz^ 



- K^B_i_ = 



(8) 



with K^ — A^^ (1 +7gp) ~ A^^. The solution reads for 
\z\ > dp/2, 



B 



_(z)=Bi(dj./2)exp. 



\z\-dF/2 



(9) 



Inserting this expression for Bj^ into Eq. (6) we obtain 
by use of Eq. (3) 
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(n^ X V±ip) ~ — Mj_ 
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where we have set Xl — Xl + dF/2. The magnetic field 
B decays exponentially with increasing z provided the 
thickness of the S layers exceeds the London penetration 
length Xl- 

In order to obtain an equation for the phase difference 
if of the superconductors we use the Maxwell equation 
(V X H)^ - c-^ d Djdt = (47r/c) j^ and the stan- 
dard expression for the Josephson current according to 
the Stewart-McCumber model^^'^^. This simple model 
[also known as the resistively and capacitively shunted 
junction (RCSJ) model] provides a good description of 
the CVC of a real J J, although effects due to finite di- 
mensions of the contacts and nonlinearities of the quasi- 
particle current are neglected. Using the Josephson rela- 
tion 



dip 



2eV 



(11) 



and the standard expression for the Josephson current 
we obtain within this model 



nf 
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- ZjV^(^ + sin((y9) = 
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V 



'2XlJ, 



(V X M^)^ (13) 



where here, too, 51 j — (2ejc/Cnfi)^/^ is the Joseph- 
son "plasma" frequency, jn — {R\jC\j)^^, Coiui) = 
e(a;)/47rd and i?n(w) are the capacitance and resistance 
of the junction per unit area, respectively, d is the thick- 



ness of the insulating layer, Ij = Vj/flj, vj — cy d/2eXL 
is the plasma wave propagation velocity (Swihart waves), 
and rj — j / jc is the normalized bias current through the 
junction. The capacitance Cn and the resistance i?n of 
the junction may depend on frequency w (in the Fourier 
representation). A simpler equation for the phase differ- 
ence ip in the stationary case has been reported previ- 
ously in Ref. 53. In a general, non-stationary case, this 
equation was derived in Ref. 30. Note that a slightly 
different approach for the study of dynamic processes in 
SFS junctions was used in a recent paper'^^. In particu- 
lar, Eq. (13) can be easily derived from Eqs. (A3)-(A6) 
of this work. 

In order to obtain a closed set of equations for the 
phase difference p of the superconductors and the mag- 
netization Mj^ of the ferromagnetic layer, we need to use 
a dynamic equation for Mj^ as well. 

The dynamics of the magnetization M in the F layer 
is described by the well-known Landau-Lifshitz-Gilbert 
(LLG) equation (see, e.g., Refs. 43 and 54), which allows 
one to describe the temporal development of M in an 
effective magnetic field Hos including all internal and 
external contributions. 

We decompose the magnetization vector M according 
to M = AfoHg -I- m, where the unit vector ng denotes 
the easy axis direction and m _L Ug is the dynamic part 
which evolves in time as described by the LLG equa- 
tion. Assuming that in equilibrium the magnetization 
coincides with the static part along the easy axis, i.e. 
Mq w |M| ^ |m|, and using B^ — 0, we obtain 



dm 
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(14) 



where M^s = g\e\/2mc, g < Qis the gyromagnetic factor, 
a is a parameter related to the anisotropy constant'*'^, Im 



is a characteristic length related to spin waves, and 7m 
is the dimensionless Gilbert damping constant. 

We further neglect the Gilbert damping term (7M = 
0), align the easy axis along the z direction (e = z), 
and substitute B_lj? = 47rMj_ + Hj_ (Mj_ = m) into 
Eq. (14), where Bj_f is the magnetic induction in the F 
layer and Hj^ is the magnetic field, which is assumed to 
be independent of the z coordinate (screening effects in 
the F layer are negligible). The field Hj_ is continuous 

across the SF interface, i.e., Hj^ = B^ ( z -^ dp 12 j , and 
is given by Eq. (10). 
Finally, we obtain 



9m 



= Vl 



M 



$0 



M X m 



ilfo 



(47r)2(a-l)A, 



-V^if 



(15) 



where VLm = 47r(a — l)|Meff|Mo is the resonance fre- 
quency of magnetic moment precession [a > \), s — 



dF/[2{a - l)Ai], ll = [al{a ~ 1)]Z|,. 



Equations (13) and (15) fully describe different dynam- 
ical processes in the junctions under consideration. Note 
that the Josephson current is coupled to the magnetiza- 
tion through the spatial derivative of the phase difference 
Vj^(^ [the last term on the right-hand side of Eq. (15)]. 
Therefore, in a spatially homogeneous case there is no 
coupling between the Josephson effect and dynamics of 
the magnetization. 



III. FISKE STEPS 

In this section, we consider a SFIFS Josephson junc- 
tion in a weak external magnetic field i^oxt assuming that 
it is constant in space and time and is directed parallel 
to the interfaces along the y direction. As is well known, 
in this case so-called Fiske steps arise on the CVC due 
to excitation of eigenmodes in the junction. The phase 
difference (^(x, t) depends on the x coordinate and, there- 
fore, dynamics of the magnetic and superfluid systems 
are coupled together. We consider the case when the 
magnetization vector in the stationary state is directed 
perpendicular to the SF interfaces, i.e. Mq — M^n^ and 
Ho = — 47rMo. As the typical values for the magnitude 
of the stationary magnetization A/q are hundreds of Gaufi 
and the small external magnetic field iJoxt is of the order 
of a few Gaufi, one can neglect the in-plane magnetization 
My — — iJcxt/(47r) compared to Mq. The resulting pre- 
cessional motion of the magnetization M in presence of a 
current through the JJ implies that the in-plane compo- 
nents m _L n^ of M are excited. Therefore, we represent 
M as M(a;,t) — Mq + m(x,i). Components nir^y are 
easily found from Eq. (15): 
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(16) 

(17) 



Equations (16), (17) are written under the assumption 
that all relevant quantities depend on time as exp(ia;i) 
and, what is more important, spatial derivatives in the 
equation for 111(2:, t) are neglected. The latter assump- 
tion is justified provided the magnetic length Im is much 
shorter than the Josephson length Ij: Im -^ h- It is 
not difficult to analyze a more general case of arbitrary 
relation between Im and Ij, but the corresponding for- 
mulas become too cumbersome. Substituting Eq. (16) 
into Eq. (13) we obtain 



uj{uj- IJr) 
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<p{x,u)+ (18) 

+T{sin{Lp)}{x,uj) = 77, 



where J^{sin((p)}(a;, w) is the Fourier transform of 
sm\ip[x^t)\ with respect to time t and 



lj{^) = h 



(1 + s)C^F 



1/2 



(19) 



is a renormalized Josephson length containing C^jF and, 
therefore, depending on frequency lo. Equation (18) is 
the favored generalization of Eq. (1) for SFIFS junctions. 
In order to find the CVC, we represent the phase 
difference tp of the superconducting layers in the form 
If = ipQ{x,t) + ij{x,t) (see Ref. 50). The first term is 
given by ipoix, t) = khx + flyt with kh = 47rALifext/^o 
[see Eq. (10)] and fly = 2eV/h. The function ip{x,t) 
is assumed to be small allowing us to linearize Eq. (13) 
with respect to ip: 



-P{ij}{x,t)=sm[^o{x,t)] (20) 

= sin(J7yi) cos{khx) + cos{Qvt) sui{khx) 

where the operator P is defined as 



p ^n, 
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The current correction drj to the dc current tjq = 

{2eV/h)/Q,j ~ VLv /^j is given by 



(577= {'ip{x,t)cos[ipo{x,t)]), 



(22) 



where the angular brackets denote the average with re- 
spect to space and time. 

Equation (22) determines the constant normaUzed 
current through the junction as a function of voltage 
V, which gives a current-voltage (I-V) curve. Equa- 
tion (20) contains parts oscillating in space and time. 
It should be solved taking into account the boundary 
conditions^-^'^^'"*^'^" 



dx 



= 0, 



(23) 



x=±L 



where L denotes the length of the junction along the x 
direction. The right-hand side of Eq. (20) can be writ- 
ten in the form Im{exp(zriy<)[cos(Kjya:;) -I- isin^KHx)]} 
and, therefore, the solution of Eq. (20) can be written as 
'ijj{x,t) = Im{exp(ir2yi)?/ji(a;)}, where the function ^'1(2;) 
obeys the equation 



- Paii'ijix) = cos{khx) + i ain{KHx) (24) 

with the boundary condition Eq. (23). 

The operator Pn coincides with P after replacing d/dt 
by iriy. The solution can be easily found and equals 



V'l — 15~7T7~E7\ 1 (^os{khx) +Ccos{kvx) +isiii{KHx) + 



+ iS sin 



{kvx)\ 



(25) 



where Pn{V,H) = a^ - l]{nv)K%, a^ = f^7^(rj^ - 
i^fiQ.v) and 
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with 9h = uhL, 6v — nyL, Ky — o?lJ {^v)- Substi- 
tuting the function ip[x,t) expressed through i^iix) into 
Eq. (22), we find the dependence, 5'q{V) = Sj{V)/jc, 



'^^'"{^^^ 
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V {ejt - ei) ' 



cos(26'v) - cos(26i/f ) 



sin(26'y) 



.(27) 



Since we assumed that the correction ■;/; = 
Imjexp (i51\/t) '0i(a;)} to the phase difference (f in 
the superconducting layers is small, Eq. (27) is only 
valid for normalized voltages fly/Jlj > (7fl/r2j)^^. 
This can be seen from Eq. (25) where one should verify 
that the prefactor P^^iV^H) is small. 

Let us discuss the current results and compare them 
with those obtained in Ref. 30. The prefactor Pq ^(V, H) 
in Eq. (27) contains the renormalized Josephson length 



I J defined in Eq. (19), which corresponds to the quantity 
Iv of Ref. 30. The formulas for Fiske steps in Ref. 30 
were given for small values of the parameter s. If the 
parameter s is not very small, one can reproduce the 
correct result by replacing there ^\j -^ (1 4- s)Vt\j, i.e., 
Eq. (19). [Note that in the definition of ^ms, Eq. (10) 
of Ref. 30, there is a misprint. The factor of two in the 
exponent at the right-hand side is missing so that the 
correct formula reads ^\[g = ^m(1 + s)"^ ] The modi- 
fied dependence of the normalized Josephson length Ij 
on the parameter s changes the form of the I-V charac- 
teristics and reveals that the effect of the ferromagnetic 
layer is much more pronounced compared to the results 
of Ref. 30 even for small s because the denominator in 
Eq. (10) of Ref. 30 is very small at voltages correspond- 
ing to peaks in the CVC and, therefore, is very sensitive 
to the parameter s. Thus, we update the figures showing 
the dependence (5ry(l^iorm.) a-s a function of normalized 
voltage T4orm. = ^v/^J- Finally, we also present I-V 
characteristics for different values of normalized junction 
lengths L/lj including those which correspond to the ex- 
perimental values of Ref. 28 {L/lj < 1). As in Ref. 30, 
for simplicity, we assume that the damping coefficient 77? 
is constant, i.e., it does not depend on voltage V. 

In Figs. 2 and 3, we plot the current correction Srj 
as a function of normalized voltage Vnorm. for different 
values of the parameter s — dF/(2(a — 1)Al) and nor- 
malized junction length L/lj. Taking into account the 
experimental values of L and Ij (see Ref. 28), we dis- 
play the current correction Sr] for short junctions with 
L/lj = 0.75 [see Figs. 2(a) and 2(b)] and, in addition, 
for longer junctions with L/lj = 2 [see Figs. 2(c) and 
2(d)] and L/lj = 10 [see Figs. 3(a) and 3(b)]. Black 
curves represent the limit s -^ where we have no F lay- 
ers in the system and the CVC correspond to ordinary 
Fiske steps. Due to the fact that in experiments, only 
the strength of the external magnetic field can be varied, 
we display our result for different values of the parameter 
khIj oc -ffoxt keeping all other system parameters such 
as flM/^j, L/lj, and s constant. 

The strongest influence of the ferromagnetic layers 
on the current- voltage characteristics develops for exter- 
nal magnetic flelds such that the parameters khIj and 
^m/^j coincide. By comparing Figs. 2(a) and 2(b) [or 
Figs. 2(c) and 2(d), respectively] one can observe that 
the change of the current correction is clearly recog- 
nizable for KhIj = ^m/^j and nonzero s, while for 
i^H^.j 7^ ^Aj/f^Ji it only becomes pronounced for larger 
values of s. 

As can be seen from Figs. 2(a) and 2(c) that the nor- 
malized junction length L/lj determines the form of the 
CVC even in the case s = 0, i.e., the number of Fiske 
steps close to the normalized magnetic resonance fre- 
quency Qm/^j may vary for different values of L/lj. 
Provided for s = there appears a single peak close to 
rim/flj, increasing the parameter s leads to a double 
splitting of the dominant peak. For even larger values 
of s, the pair of peaks moves more and more apart from 
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FIG. 2. (Color online) Correction to the I-V characteristics of an SFIFS JJ in a weak external magnetic field due to interaction 
of Josephson oscillations and spin- wave modes. The correction is plotted as a function of normalized voltage Kiorm. = fiv/^j 
for different values of the parameter s. The figures are presented for the following parameters: (a) Q,m l^j = i^h^j = 8 and 
L/lj = 0.75, (b) nAi/^j = 8, K„lj = 12, L/lj = 0.75; (c) Qm/^j = k.hIj = 8, L/lj = 2; (d) flAi/^j = 8, k„Ij = 12, L/lj = 2. 
The damping coefficients are "/r/Qj = 0.4, 7A//nj = 0.3. 



each other [see Fig. 2(a)]. A similar effect can be seen 
for a larger number of Fiske steps close to fin/j/^j, e.g., 
Fig. 2(c) displays essentially two Fiske steps in the vicin- 
ity of VIm/^j = 8 that both split up into two peaks 
moving apart from each other with increasing s. 

For distinct values of the parameters khIj and nm/Clj 
[see Figs. 2(b) and 2(d)], there also emerge additional 
peaks in the I-V characteristics close to the normalized 
magnetic resonance frequency, but the detailed impact of 
the F layers on the CVC is not as obvious as is the case 
for Knh = ^m/^j- From Fig. 2(d), one can already con- 
jecture that for long junctions, the ferromagnetic layers 
simply induce a single additional peak close to VIm/^j- 
In Fig. 3, where the current correction Srj is shown for 
the limit of large values of L/lj {L/lj = 10), this fea- 
ture becomes more apparent. For coinciding values of the 
magnetic resonance frequency JIm/^j and the parame- 
ter Kfjlj [see Fig. 3(a)], we find a single peak for s = 
and a double peak for s 7^ in the vicinity of flM/^.i- 
For fluf/^j y^ khIj, there emerges a single peak close to 
^Im/^j and krIj, respectively, where the former is no- 
tably smaller in magnitude [see Fig. 3(b)]. Below we also 
derive analytical expressions for these peak positions. 

Thus the presence of the F layers leads not only to a 



shift of the peaks in the dependence 5?7(V^oim.) but also 
to a change of the overall form of this dependence. The 
additional peaks arising on the I-V curves can be at- 
tributed to the ferromagnetic resonance and the nonzero 
coupling between Josephson and magnetic moment oscil- 
lations. In order to observe these peaks experimentally, 
one should perform measurements with different samples 
that contain ferromagnetic layers of varying thickness. 
Then, according to our theoretical result, one would be 
able to differentiate between ordinary Fiske steps and 
peaks caused by interaction of Josephson and magnetic 
oscillations in the F layers. Note that in the limit of a 
very short junction {L/lj <?C 1) there is no coupling be- 
tween Josephson and magnetic moment oscillations. In- 
deed, in this limit we obtain from Eq. (27) 



Sr] ~ Im 



^j 



fly i^V - Hr) 



(28) 



It is seen that magnetic characteristics such as J7j\/ of the 
F layers drop out from this expression. 

In the limit of long junctions, L/lj ^ 1, the expression 
for the current correction can be approximated by 
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FIG. 3. (Color online) Current correction Srj as a function of 
normalized voltage Vnorm. ~ Q,v/^J for long JJ with normal- 
ized junction length L/lj = 10. The function Sr] is displayed 
for (a) Qm/^j = i^hIj = 8; (b) Qm/Oj = S,k„Ij = 12 
and for different values of the parameter s. The damping 
coefficients are "fn/i^j = 0.4, 7A//r2,/ — 0.3. 
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Pn{V,H) 



(29) 
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In accordance to Fig. 3 we obtain for s = a single peak 
at normalized voltage Vnorm. — i^h^j while for s 7^ and 
khIj = rij\//r2j there exist two peaks at 



K, 



VL 



M 



n,, 



^{i + s) ± ys(TT^. 



(30) 



Finally, for the general case s 7^ and khIj 7^ flM/^j 
we find in leading order in the parameter s two peaks 
located at normalized voltages 



^norm. 




1 x^ 

/ 1 1 " 


(31a) 


r ' 1-X2 
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(31b) 


•^norm. 


i 1 J ' „ 
1 — x"^ 



where a; = {VLm /^j)/{k.hIj). These analytical expres- 
sions perfectly describe the peak locations of the current- 
voltage characteristics in the limit L/lj ^ 1 as exemplar- 
ily shown in Fig. 3 for junctions with L/lj = 10. 



IV. COUPLED COLLECTIVE MODES 

In this section, we analyze the spectrum of coupled 
collective modes in long Josephson junctions with a fer- 
romagnetic layer. So far we have derived essentially two 
(coupled) equations, Eqs. (13) and (15), that describe re- 
spectively the dynamics of the phase difference </? of the 
S layers and the magnetization M of the ferromagnetic 
layers. Here, we consider again the case when the mag- 
netization Mo is aligned normal to the interface so that 
in equilibrium Bq — 0. Small perturbations near the 
equilibrium result in precessional motion of the magnetic 
moment M and in a variation of the phase difference 
Lp in space and time. In order to find the spectrum of 
collective modes in the system, we represent the phase 
difference (p and the magnetic moment M in the form 



<p = (^o + '0, M = Afon^ -Fm_L, 



(32) 



where n^ is the unit vector normal to the SF interface 
and the functions ^ and m_L are assumed to be small, 
IV^I ^ Wo\ and |m_L| <C |Mo|. Linearizing Eq. (13) with 
respect to "0, we find that the function ip{x,t) obeys the 
equation 



n-^ 



2 f(^ 



IR 



dt 



-fjWU- 



-V 



(33) 



cdf 



2XlJc 



V X m_L 



The perturbation mj^ of the magnetic moment is parallel 
to the SF interface and is described by the equation 



dni± 
dt 



= n 




lliVl) 



n^ X mj^ 



(34) 



$n 



(4^)2/3A. 



-ViV^ 



-7M 



dt 



where we included again the Gilbert damping term, 
which was neglected in Eq. (15). Fourier transforming 
the perturbations (p{r,t) and m_L(r,t) to (k_L,w) repre- 
sentation and combining Eqs. (33) and (34) into a single 
equation, we obtain 



Ml ^^^^'"\ 1=0, 



'n 



M 



J (wj — w^) ibky —ibkx 



ZCCKy 



-COM 






(35) 
(36) 
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In order to show this exphcitly, we consider the case 
without damping, 77^ = 7m = 0, and assume that s <C 1 
and Im ^ ^j, which means that we neglect the spatial 
dispersion of spin waves on the Josephson length (these 
conditions are usually fulfilled experimentally) . It is con- 
venient to write Eq. (37) in the dimensionless form 



[Z - 1 - g2] \Z - Zm] = sq^Z, 



Ml 



(38) 



where Z = (uj/Qj)^, Zm = {^m/^jY ^-i^d q = Ijk. One 
can see that for s = the two dispersion curves Z = l+q"^ 
and Z = Zm cross each other at q^ = Zm — 1- To find 
the form of the dispersion curve in the vicinity of the 
crossing point q^, we represent Z and q, respectively, as 
Z = Zm + 5Z and q — qo + Sq. Then, one can easily 
obtain from Eq. (38) 



SZ = qo 



Sq± y/Sq 



sZ 



M 



(39) 



In Fig. 4 we plot the spectrum of coupled spin and 
plasma-like modes and the function SZ(Sq) close to the 
crossing point. Here, we take into account a finite value 
of the parameter Im/Ij so that Eq. (38) that determines 
the function Z{q) takes the form 



FIG. 4. (Color online) Spectrum of coupled spin and plasma- 
like modes, (a) The function Z{q) = {uj{q)/ilj)'^ determined 
by Eq. (40) is shown for finite values of s and in the inset for 
the case s = 0. (b) The dependence 5Z{5q) that represents 
the spectrum close to the crossing point is plotted according 
to Eq. (39). The following parameters are chosen: Zm = 
4:, Im/Ij ^0.1. 



{Q,M/njf 



,j = u!'j{k,u!) = ilj (l + fc^/j) - i'juuj, um = 
ujMik,uj) = Q,M (1 + s + k'^llj) - i-iM^i b = sj3c/jc, a = 
sflMl^/b, and k — |kj^|. 

The homogeneous equation (35) has a non-vanishing 
solution provided the determinant of M equals zero. Set- 
ting det(A^) equal to zero we obtain the dispersion rela- 
tion 



[z- 



z-z 



M 



sq 



ZmZm, 



(40) 



with Zm = ZM[l+sW{hdlljfY- The inset of Fig. 4(a) 
indicates that the two branches indeed cross each other 
for s = 0, whereas for s 7^ we find a "repulsion" of the 
spin and Josephson excitations. Figure 4(b) displays the 
function 5Z{5q) that represents the behavior of the spec- 
trum in the vicinity of the crossing point and distinctly 
emphasizes the mutual repulsion. Both the dispersion 
curves Z{q) and 5Z{5q) given by Eqs. (39) and (40), re- 
spectively, are presented for several values of s and the 
parameters JIm/^j = 2, Im/Ij = 0.1. 



FERROMAGNETIC RESONANCE 



;j] [w^ - ujm] = svjilMUJMk'^ 



(37) 



From Eq. (37) we can conclude that the spin and 
charge excitations decouple only in the limit when the 
right-hand side of this equation can be neglected. In 
this case the spin waves with spectrum wa/ (fc, w) and the 
plasmalike Josephson waves with spectrum ujj{k,ijj) ex- 
ist separately. In the general case, Eq. (37) describes the 
spectrum of coupled spin waves and plasma-like modes in 
the system. The most interesting behavior corresponds 
to the case flM > ^j- In this situation, the two branches 
of the spectrum cross each other in the absence of the cou- 
pling, while a finite coupling leads to mutual repulsion of 
these branches. 



In this section, we study the response of the sys- 
tem to an external oscillating magnetic field 7fext(0 = 
Hi, sm{iyt) with a small amplitude Hi, ^ Mq and fre- 
quency I'. The applied field is supposed to be directed 
along the y axis, i.e., Hoxt(0 = -ffoxt(0'^y We as- 
sume again that the equilibrium magnetization Mq is 
oriented in the z direction, M — Mo^^z- The external 
magnetic field -ffcxt (t) causes precessional motion of the 
magnetization vector M and a variation of the phase dif- 
ference ip in space and time. As before (see Sec. Ill), 
we, respectively, represent magnetization and phase dif- 
ference in the form M(a;,t) = MqUz + ■my{x,t)ny and 
ip{x, t) = ipQ + 'ipiXi t). Here, (po is a constant determined 
by a bias current jb = jc sin((/5o) and ipi^i t)i iTT-yi^, t) are 
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FIG. 5. (Color online) Frequency dependence of the normalized absorbed power Vn in the JJ as a function of the normalized 
frequency uj/Q,j for different values of the parameter s and for different damping coefficients jr/Qj and jm/^j- The figures (a) 
and (b) are shown for the parameters L„ — 3, yM/^j = 0.1, (fiM/fij)^ — 5 and, respectively, (a) jn/^j = 0.1, (b) 'Jr/Qj — 1. 
With regard to figures (c) and (d) we have chosen Ln = 1, {0,i\i/fl,j)'^ = 5 and in (c) jm/^j ~ 0.1. 



small perturbations due to the external ac magnetic field 
Ho^tit) = H^lm[exp{iiyt)], IV'I <C \^o\, \my\ < |Mo|. 
Due to the coupling of M and (p, we expect modifications 
of the ferromagnetic resonance in the system appearing 
as additional features in absorption spectra. 

Thus, to study ferromagnetic resonance, we need to 
calculate the power V (per unit area) absorbed in the 
system. The absorbed power V can be found as the time- 
averaged difference between the energy flux 5in,out com- 
ing in and out of the system. These fluxes are expressed 
in terms of Poynting vectors S^^ 



applied weak ac magnetic field iJcxt (t) ■ This vector S dif- 
fers from zero only in the insulating layer of thickness d. 
The magnetic field consists only of the applied ac field, 
H(a; — ±L) = i?ext(i)ny and, therefore, the Poynting 
vectors are directed parallel to the x— axis. We repre- 
sent the phase difference in form of the Fourier trans- 
form ip(x,t) — J du!/{2Tr)e:xjp{iut)'ip{x,Uj), iJj{x,—U!) = 
iIj*{x,uj). 

The function tp^XjCo) obeys an equation that is derived 
in a way similar to the derivation of Eqs. (18)-(22) and 
has the form 



"P = / dzdy n^ • (5i„ - So 



(41) 



where 5in,out = (c/47r)[E x H]2;=±l and the angular 
brackets denote averaging with respect to time t. 

The electric field E — UzE is directed along the z— axis 
and is related to the time derivative of the phase differ- 
ence via the Josephson relation 



E = -{l/d){h/2e)d^/dt. 



(42) 



Therefore, in order to find the Poynting vector 5, we have 
to calculate the function ipit) which is determined by an 



dxl 



1p{Xnj'jj) - Kl-ipiXn,U;) = 0, 



(43) 



where we have introduced the dimensionless variable 
Xn = x/lj and have set 



£. 



^ujj 



1 + .s/[(l + s)C^f] 



(44) 



withCi^j = cos{ipo) ~ ijj{Lu-i^ Ft)/ ^"^j, a^j = l+s/[(l + 
s)Ci^f], and £„f is defined in Eq. (17). Equation (44) is 
supplemented by the boundary conditions 
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FIG. 6. (Color online) Frequency dependence of the normalized absorbed power Vn as a function of the normalized frequency 
uj/Q.j for really short junctions. Here, too, the figures are displayed for different values of the parameter s and for the choice 
^r/^j ~ l,7A//f2./ = 0.1, {Q,m/Q,j)'^ — 5. In figure (a) L„ — 0.2 whereas in figure (b) L„ = 0.6. 



a-ui-^ — -0(2;, w) 



x=±L 



and Hq = — ^ 



(45) 
that can be obtained from Eqs. (9) and (10). As a con- 
sequence, the solution for Eq. (43) has the form 



n±L, u) = T Tt n (46) 

where O^j — K^i„ = 9'^ + i9'^, Ln = L/lj. Fourier trans- 
forming Eq. (46) back into the time representation we 
obtain 



no 



tanh(6'!y) 



(47) 



Taking into account that all the quantities do not depend 
on y, the absorbed power V can be represented as 



he l di,[L,t) 



(48) 



where Ly is the length of the junction along the y direc- 
tion. Substituting Eq. (47) into Eq. (48) and relabeling 
the external field frequency i^ — )■ w, we finally arrive at 



V = 



$0^0 (n. 



(2^)2 \Ho 



LxLy 

— ; — -uj Im 



tanh {9i^) 



(49) 



where L^ = L. This formula differs drastically from a 
standard formula for the absorbed power V in ferromag- 
netic films because it describes the power absorption not 
only in the F film, but also in the Josephson junction. In 
particular, V y^ even in the absence of the ferromag- 
netic layer. In this case, Eq. (49) describes the power 
needed to excite standing plasma waves. 



In Fig. 5, we plot the frequency depen- 
dence of the normalized absorbed power Vn = 
(a;/il,/)Im [tan]i{9i^) / (ai^9^)] as a function of normalized 
frequency w/fij at different s and normalized junction 
length L/lj. Generally speaking, from Fig. 5(a), we see 
that at s = (no ferromagnetic layer), there are periodic 
resonances related to excitation of standing waves in 
the Josephson junction (Josephson plasma resonances). 
Interestingly, in the presence of the ferromagnetic layers, 
s 7^ 0, additional peaks appear on the curves. These 
peaks are caused by the ferromagnetic resonance in 
the F layer at frequencies uj sa D,m- With increasing 
s the influence of the F layer becomes more and more 
pronounced. One can see this from the fact that, for 
instance, the spectrum close to a; « ^Im appears to have 
a more complicated structure and the peaks increase in 
height. 

To indicate the influence of the (normalized) damp- 
ing parameters 7^ and 7m, we display in Fig. 5(b) the 
normalized absorbed power P„ for the case 7/? >• 7^/. 
This shows that the periodic resonances in the junction 
are strongly suppressed and the absorption spectrum is 
dominated by the effect of the ferromagnetic layer. In 
addition to that, the normalized length of the junction 
Ln = Lx/lj also determines the absorption spectrum. In 
Figs. 5(c) and 5(d), Pn is shown for the case L„ = 1 
and here, too, for different values of the parameters s,^b., 
and 7m- We find that the distance between periodic res- 
onances is larger for short junctions and compared to 
Fig. 5(a), where L„ = 3, the infiuence of the F layer 
on the absorption spectrum is weaker. The blue curve 
in Fig. 5(c) reveals that the periodic resonances can be 
almost completely suppressed by increasing the damping 
parameter 7^'. Eventually, Fig. 5(d) indicates that in 
systems where both damping parameters ^n and ji\i are 
large and of the same order of magnitude, the effect of 
the ferromagnetic layer becomes negligible. 

In Fig. 6, we show the frequency dependence of the ab- 
sorbed power Pn for short junctions of length L„ = 0.2 
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[see Fig. 6(a)] and L„ = 0.6 [see Fig. 6(b)]. In Fig. 6(a), 
the peak at co/flj « 2.3 is related to the ferromagnetic 
resonance in the F film. In contrast to this, slightly longer 
junctions feature a much stronger influence of the fer- 
romagnetic layer as becomes apparent from Fig. 6(b). 
More importantly, we find that the relative magnitudes 
of peaks due to the Josephson plasma resonances and the 
ferromagnetic resonance are even in the case of 7/j 3> 7_a/ 
of the same order of magnitude for short junctions con- 
trary to Fig. 5(b), L„ = 3, where the Josephson plasma 
resonances are considerably smaller for the same choice 
of parameters. 

Note that our analysis is valid for not too high fre- 
quencies as it is assumed that the penetration depth is 
not frequency dependent. This means that the inequal- 
ity uj <^ ^jIj/Xl — vj/\l should be fulfilled. For this 
reason. Figs. 6(a) and 6(b) indicate only a small number 
of Josephson plasma resonances. 



VI. DISCUSSION 

We studied dynamic properties of Josephson junctions 
with a magnetically active layer characterized by the 
magnetic susceptibility x(li;,A:). These junctions may be 
of the SFIFS or SIFS type with conducting or insulat- 
ing ferromagnets. In the former case, we assumed that 
both vectors Mi and M2 characterizing the stationary 
orientation of magnetization in the F layers were aligned 
along the z direction, and our results are applicable only 
in this situation. 

We calculated the form of the CVC for SFIFS junctions 
in the presence of a weak magnetic field and found a 
modification of Fiske steps due to the presence of the 
ferromagnetic layer. The position of these steps depends 
on the relation between different parameters, especially 
between khIj and ^Im/^j- 

We have also analyzed the spectrum of the collective 



coupled modes in long JJs with a ferromagnetic layer. 
If the frequency of the ferromagnetic resonance 51m is 
higher than the characteristic Josephson frequency flj, 
then coupled magneto-plasma modes (spin waves and 
Josephson plasma-like modes) occur in the region of 
crossing terms. 

The analysis of the ferromagnetic resonance in the F 
layer incorporated in JJs of the SFS or SFIFS types shows 
that the peaks in the frequency dependence of the ab- 
sorbed power P{uj) correspond both to the ferromagnetic 
resonance in the F film and to the Josephson plasma res- 
onances in the tunnel J J. 

It is not easy to compare our results with available 
experimental data. The dynamic properties of ferromag- 
netic layers play a crucial role in determining the form of 
the CVC (Fiske steps). Meanwhile, little is known about 
these properties in experiments. It would be useful to 
study experimentally magnetic resonance in the F lay- 
ers at temperatures above the critical temperature of the 
superconducting transition Tc- The frequencies of the 
Josephson oscillations Qj and magnetic resonance ^m 
should not be very different. In addition, we assumed 
that the easy-axis magnetization is perpendicular to the 
SF interface. There are no data about magnetization 
orientation in junctions studied experimentally. 

As to magnetic resonance, we are only aware of 
Refs. 32, 55, and 56 where ferromagnetic resonance was 
measured on SF structures. However, the authors of 
Ref. 32 measured the CVC of a SFS junction with a 
strong damping, but not the absorbed power. In Ref. 55 
and 56, the absorbed power was measured, however not 
in SIFS junctions, but in SF bilayers. Thus further 
experiments are needed to study the interplay between 
magnetic and Josephson oscillations in tunnel Josephson 
junctions with a ferromagnetic layer. 

ACKNOWLEDGMENT 

We thank SFB 491 for financial support. 



^ A. A. Golubov, M. Y. Kupriyanov, and E. Il'ichev, Rev. 

Mod. Phys. 76, 411 (2004). 
2 A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 
^ F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. 

Phys. 77, 1321 (2005). 

* M. Eschrig, Physics Today 64, 43 (2011). 

^ L. N. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, Sov. 

Phys. JETP Lett. 25, 290 (1977). 
® A. I. Buzdin, L. N. Bulaevskii, and S. V. Panyukov, Sov. 

Phys. JETP Lett. 35, 178 (1982). 
'' V. V. Ryazanov, V. A. Oboznov, A. Y. Rusanov, A. V. 

Veretennikov, A. A. Golubov, and J. Aarts, Phys. Rev. 

Lett. 86, 2427 (2001). 

* V. A. Oboznov, V. V. Bol'ginov, A. K. Feofanov, V. V. 
Ryazanov, and A. I. Buzdin, Phys. Rev. Lett. 96, 197003 
(2006). 

® T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stephanidis, 



and R. Boursier, Phys. Rev. Lett. 89, 137007 (2002). 
^° Y. Blum, A. Tsukernik, M. Karpovski, and A. Palevski, 

Phys. Rev. Lett. 89, 187004 (2002). 
^^ A. Bauer, J. Bentner, M. Aprili, M. L. DeUa Rocca, 

M. Reinwald, W. Wegscheider, and C. Strunk, Phys. Rev. 

Lett. 92, 217001 (2004). 
^^ H. Sellier, C. Baraduc, F. Lefloch, and R. Calemczuk, 

Phys. Rev. Lett. 92, 257005 (2004). 
^^ V. Shelukhin, A. Tsukernik, M. Karpovski, Y. Blum, 

K. B. Efetov, A. F. Volkov, T. Champel, M. Eschrig, 

T. Lofwander, G. Schon, et al., Phys. Rev. B 73, 174506 

(2006). 
^'^ A. A. Bannykh, J. Pfeiffer, V. S. Stolyarov, I. E. Batov, 

V. V. Ryazanov, and M. Weides, Phys. Rev. B 79, 054501 

(2009). 
^^ M. Eschrig, T. Lofwander, T. Champel, J. Cuevas, 

J. Kopu, and G. Schon, J. Low Temp. Phys. 147, 457 



13 



(2007). 
^^ A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 

(2003). 
^■^ R. S. Keizer, S. T. B. Goennenwein, T. M. Klapwijk, 

G. Miao, G. Xiao, and A. Gupta, Nature 439, 825 (2006). 
^® I. Sosnin, H. Cho, V. T. Petrashov, and A. F. Volkov, Pliys. 

Rev. Lett. 96, 157002 (2006). 
1^ T. S. Khaire, M. A. Khasawneh, W. P. Pratt, and N. O. 

Birge, Phys. Rev. Lett. 104, 137002 (2010). 
2° J. Wang, M. Singh, M. Tian, N. Kumar, B. Liu, C. Shi, 

J. K. Jain, N. Samarth, T. E. Mallouk, and M. H. W. 

Ghan, Nature Physics 6, 389 (2010). 
^^ J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, 

Science 329, 59 (2010). 
^^ D. Sprungmann, K. Westerholt, H. Zabel, M. Weides, and 

H. Kohlstedt, Phys. Rev. B 82, 060505 (2010). 
^^ M. S. Anwar, F. Czeschka, M. Hesselberth, M. Porcu, and 

J. Aarts, Phys. Rev. B 82, 100501 (2010). 
^^ F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 

Lett. 86, 3140 (2001). 
'^^ J. W. A. Robinson, G. B. Halasz, A. I. Buzdin, and M. G. 

Blamire, Phys. Rev. Lett. 104, 207001 (2010). 
^® M. Weides, M. Kemmler, H. Kohlstedt, R. Waser, 

D. Koelle, R. Kleiner, and E. Goldobin, Phys. Rev. Lett. 
97, 247001 (2006). 

^"^ J. Pfeiffer, M. Kemmler, D. Koelle, R. Kleiner, 

E. Goldobin, M. Weides, A. K. Feofanov, J. Lisenfeld, and 
A. V. Ustinov, Phys. Rev. B 77, 214506 (2008). 

^* Wild, G., Probst, C., Marx, A., and Gross, R., Eur. Phys. 

J. B 78, 509 (2010). 
2^ M. Kemmler, M. Weides, M. Weiler, M. Opel, S. T. B. 

Goennenwein, A. S. Vasenko, A. A. Golubov, H. Kohlstedt, 

D. Koelle, R. Kleiner, et al., Phys. Rev. B 81, 054522 

(2010). 
^° A. F. Volkov and K. B. Efetov, Phys. Rev. Lett. 103, 

037003 (2009). 
^^ S. ichi Hikino, M. Mori, S. Takahashi, and S. Maekawa, J. 

Phys. Soc. Jpn. 77, 053707 (2008). 
^^ L Petkovic, M. Aprih, S. E. Barnes, F. Beuneu, and 

S. Maekawa, Phys. Rev. B 80, 220502 (2009). 
^^ X. Waintal and P. W. Brouwer, Phys. Rev. B 65, 054407 

(2002). 
^^ V. Braude and Y. M. Blanter, Phys. Rev. Lett. 100, 

207001 (2008). 



^^ E. Zhao and J. A. Sauls, Phys. Rev. B 78, 174511 (2008). 
'^^ F. Konschelle and A. Buzdin, Phys. Rev. Lett. 102, 017001 

(2009). 
^^ J.-X. Zhu, Z. Nussinov, A. Shnirman, and A. V. Balatsky, 

Phys. Rev. Lett. 92, 107001 (2004). 
^* J. Michelsen, V. S. Shumeiko, and G. Wendin, Phys. Rev. 

B 77, 184506 (2008). 
^^ G. Holmqvist, S. Teber, and M. Fogelstrom, Phys. Rev. B 

83, 104521 (2011). 
""^ I. O. Kulik and I. K. Janson, The Josephson Effect in 

Superconductive Tunneling Structures (Israel Program for 

Scientific Translations, Jerusalem, 1972). 
*^ K. K. Likharev, Dynamics of Josephson Junctions and 

Circuits (Gordon and Breach Science Publishers, London, 

1986). 
"^^ A. Barone and G. Paterno, Physics and Applications of the 

Josephson Effect (John Wiley & Sons, New York, 1982). 
*^ L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Elec- 
trodynamics of Continuous Media, Course of Theoretical 

Physics Vol. 8 (Butterworth-Heinemann, Oxford, 1984). 
■'■* B. D. Josephson, Rev. Mod. Phys. 36, 216 (1964). 
■'^ B. D. Josephson, Advances in Physics 14, 419 (1965). 
■"^ A. Buzdin and I. Baladie, Phys. Rev. B 67, 184519 (2003). 
^■^ A. S. Vasenko, S. Kawabata, A. A. Golubov, M. Y. 

Kupriyanov, C. Lacroix, F. S. Bergeret, and F. W. J. 

Hekking, Phys. Rev. B 84, 024524 (2011). 
■'^ N. G. Pugach, M. Yu. Kupriyanov, E. Goldobin, 

R. Kleiner, and D. Koelle, Phys. Rev. B 84, 144513 (2011). 
^^ I. O. Kulik, Sov. Phys. JETP Lett. 2, 84 (1965). 
^° R. E. Eck, D. J. Scalapino, and B. N. Taylor, Phys. Rev. 

Lett. 13, 15 (1964). 
^^ W. C. Stewart, Appl. Phys. Lett. 12, 277 (1968). 
'^^ D. E. McCumber, J. Appl. Phys. 39, 3113 (1968). 
^^ A. F. Volkov and A. Anishchanka, Phys. Rev. B 71, 024501 

(2005). 
^* A. Aharoni, Introduction to the Theory of Ferromagnetism, 

International Series of Monographs on Physics (Oxford, 

Glarendon Press, 1996). 
^^ I. GarifuUin, D. Tikhonov, N. Garif'yanov, M. Fattakhov, 

K. Theis-Brohl, K. Westerholt, and H. Zabel, Apphed 

Magnetic Resonance 22, 439 (2002). 
'^^ G. Bell, S. Milikisyants, M. Huber, and J. Aarts, Phys. 

Rev. Lett. 100, 047002 (2008). 



